Unifying model of driven polymer translocation 
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We present a Brownian dynamics model of driven polymer translocation, in which non-equilibrium 
memory effects arising from tension propagation (TP) along the cis side subchain are incorporated 
as a time-dependent friction. To solve the effective friction, we develop a finite chain length TP 
formalism, based on the idea suggested by Sakaue [Sakaue, PRE 76, 021803 (2007)]. We validate 
the model by numerical comparisons with high-accuracy molecular dynamics simulations, showing 
excellent agreement in a wide range of parameters. Our results show that the dynamics of driven 
translocation is dominated by the non-equilibrium TP along the cis side subchain. Furthermore, by 
solving the model for chain lengths up to 10^° monomers, we show that the chain lengths probed 
by experiments and simulations are typically orders of magnitude below the asymptotic limit. This 
explains both the considerable scatter in the observed scaling of translocation time w.r.t. chain 
length, and some of the shortcomings of present theories. Our study shows that for a quantitative 
theory of polymer translocation, explicit consideration of finite chain length effects is required. 

PACS numbers: 



I. INTRODUCTION 



The transport of a polymer across a nanoporc is vi- 
tal to many biological processes, such as DNA and RNA 
translocation through nuclear pores, protein transport 
across membrane channels and virus injection Due to 
various potential technological applications such as rapid 
DNA sequencing, gene therapy and controlled drug de- 
livery 0, polymer translocation has received consider- 
able experimental 0-I3l and theoretical interest [5|-|28|. 
Of particular technological interest is the case of driven 
translocation, where the process is facilitated by an ex- 
ternal driving force. The key theoretical issue here is to 
find a unifying physical description that yields the correct 
dynamical behavior, e.g., the dependence of the translo- 
cation time T on the chain length Nq. Experiments and 
numerical simulations have indicated that r oc Nq . How- 
ever, numerous different values of a have been observed, 
suggesting explicit dependence on the various physical 
parameters (cf. Ref. Q for a recent review) . Several the- 
ories of driven polymer translocation have emerged 
[l6j . some claiming agreement with the experimental or 
numerical results within a certain subset of the physical 
parameter space. However, to date no single theory has 
been able to capture the wide range of observed values of 
a, nor quantitatively explain the reason for their depen- 
dence on the system's parameters. Therefore, the need 
for a unifying theory of driven translocation remains. 

In Refs. 0, 0], polymer translocation was described 
as a one-dimensional barrier crossing problem of the 
translocation coordinate s (the length of the subchain 
on the trans side). Here, the chain starts from the cis 
side with one end inside the pore (s = 0) and is con- 
sidered as translocated once s = uNq, with a the seg- 



ment length. The free-energy due to chain entropy and 
the chemical potential difference A/i is J^{s) = (1 — 
7')fcBTln [f {No - f )] + f A^. Here 7' is the surface ex- 
ponent (7' = 0.5, ss 0.69, « 0.95 for an ideal chain, self- 
avoiding chain in 2D and 3D, respectively), and fc^T is 
the thermal energy. From J^(s), the Brownian dynamics 
equation for s in the overdamped limit follows as = 



(l-7')fcBT 



2N0-S 



^ -hC(i)- Here F is the (con- 



stant) effective friction, and (^(t) is Gaussian white noise 
satisfying (C(t)) = and {C{t)C{t')) = 2TkBT5{t~t'). For 
moderate to large A^, this model describes translocation 
at constant mean velocity {ds/dt) = — A/i/aF. However, 
it is known that the translocation process initially slows 
down and finally speeds up towards the end 0, [3412] ■ 
Qualitatively, this observation has been explained by a 
simple force-balance argument, where the friction F de- 
pends on the number of moving monomers on the cis 
side subchain % [TUll. In Refs. [1, [H [H, it is as- 
sumed that the whole subchain on the cis side is set into 
motion immediately after the force at the pore is applied. 
However, this assumption is only valid in the limit of ex- 
tremely weak driving force, where the cis side subchain 
is always at equilibrium. In most cases, the driving force 
is substantially larger, im plyi ng that even the subchains 
are out of equilibrium [2ll424| . It has been proposed that 
in this regime, the out-of-equilibrium dynamics can be 
described by tension propagation (TP) along the chain 
backbone, which leads to nontrivial time-dependence of 
the drag force and gives the non-monotonic translocation 
velocity (TT|-[T5|. However, this idea has not been quan- 
titatively verified, since most of the studies have consid- 
ered the asymptotic limit Nq ^ 00, which is out of reach 
of experiments and numerical simulations. Therefore, it 
is imperative to study the TP mechanism for finite A'o, 
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which is the regime that is experimentally relevant and 
where numerical simulation data are available. 

To this end, in this work we adopt the TP formal- 
ism in context of the Brownian dynamics (BD) equa- 
tion of motion for s mentioned above, in which wc in- 
troduce a time-dependent friction coefficient F ~ r(t) 
that is determined by the TP equations. We introduce 
a TP formalism for finite chain lengths by incorporating 
the pore-polymer interactions to the TP equations. We 
solve the resulting Brownian dynamics - tension prop- 
agation (BDTP) model at finite chain length iVo, and 
validate it through extensive comparisons with molec- 
ular dynamics (MD) simulations. We verify that the 
tension propagation mechanism dominates the dynam- 
ics of driven translocation. In addition, we show that the 
model quantitatively reproduces the numerical values of 
a in various regimes without any free parameters, ex- 
plaining the diversity in a as a finite chain length effect. 
Finally, we address the recent theoretical disagreement 
between the constant- velocity TP theory of Refs. [TTI - [l3 | 
and the constant-flux TP theory of Rcf. [l5| and show 
that at the asymptotic limit, A^o ^ cxo, a approaches 
a = \ ^ V . 



II. MODEL 

A. General formulation 

We introduce dimensionless units for length, force, 
time, velocity and friction as s = s/a, / = /a/fc^T, 
t = iksT /rjo^ , v = vija/ksT and F = F/77, where 77 is 
the solvent friction per monomer. In these units, the BD 
equation reads 



ft>0 



f(«)§ = (i-y) 



1 



No 



+ / + CW^/tot, (1) 



where, for simplicity, we have assumed that the pore 
length Ip = a. Generalization of Eq. ([T]) to different pore 
lengths is straightforward (see, e.g., Ref [l^)- Eq. ([!]) is, 
of course, approximative rather than rigorously exact. It 
contains two approximations, which wc will show to be 
valid by quantitative comparison with MD simulations. 
First, we postulate that the friction T{t) is determined by 
TP on the cis side subchain. While there is no conclusive 
a priori reason to neglect the non-equilibrium effects of 
the trans side subchain, we will show that those effects are 
negligibly small in the experimentally and computation- 
ally relevant regimes. Second, we note that Eq. ([T]) in- 
cludes the entropic force, whose form is strictly valid only 
for small driving forces /, when the translocation time f 
is comparable to the Rouse relaxation time [H, HO] ■ 
However, for larger /, the average contribution of the 
entropic force to the total force /tot is very small (see 
results below for discussion). Therefore, even for large 
forces, the model will be shown to give excellent agree- 
ment with MD simulations. 




FIG. 1: (Color online) A 
snapshot of a translocating 
polymer in a stem-flower con- 
flguration. A tension front 
at i = —R (black arc) di- 
vides the chain into moving 
and nonmoving parts, with 
the last moving monomer la- 
beled as N. The number of 
translocated monomers is s. 



The effective friction f (f) actually consists of two con- 
tributions. The first one is the drag force of the cis side 
subchain that is solved with the TP formalism. The other 
one is the frictional interaction between the pore and the 
polymer. Formally, we can write F as the sum of the 
cis side subchain and pore frictions, T{t) — ?ycis(^) + Vp- 
While for Nq —> 00 the first term dominates, for finite 
A''o the pore friction can significantly affect the translo- 
cation dynamics. We will come back to this issue later, 
but let us first look at how the time-dependent part of 
the friction can be determined from the TP formalism. 
In the special case of extremely large driving force, one 
can find r^cis directly from the TP equations. More gen- 
erally, however, it is easier to derive the velocity of the 
monomers at the pore entrance, vq. In such a case, the 
effective friction is in a natural way defined as 



^o(i)wo(?)' 



(2) 



where ctq is the line density of monomers near the pore 
and ctqWo = ds/dt is the fiux of monomers through the 
pore entrance. In either case, determining T{t) essentially 
reduces to calculating the number of moving monomers, 
whose combined drag force then constitutes the time- 
dependent part of the friction. As the driving force is 
applied, the chain begins to move in stages, with the 
segments closest to the pore being set into motion first. 
A close analogue is a coil of rope pulled from one end, 
which first uncoils before starting to move as a whole. 
To keep track of the moving part of the chain, one de- 
fines a tension front, which divides the chain into the 
moving part that is under tension, and the nonmoving 
part outside the front (sec Fig. [1]). The front is located 
at i = —R(t), and propagates in time as parts of the 
chain further away from the pore are set in motion. The 
last monomer within the tension front is labeled as N{t). 
Using the TP formalism, one can derive an equation of 
motion for the tension front, using either R 01 N as the 
dynamical variable. The details of this calculation can 
be found in the Appendix. 



B. Different regimes 

Depending on the magnitude of the driving force, the 
equation of motion for the tension front attains a slightly 
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different form. In the simplest case, when the driving 
force is very large compared to temperature and chain 
length, / > iVg . the moving part of the chain is almost 
completely straight. In this strong stretching (SS) regime, 
the equation of motion is [lH 



dt 



fu 



T(t) (1 - z/A,7V--i) 



(3) 



Here, T{i) = N{i) - + with N{i) - S{i) being 
the number of moving monomers on the cis side. The 
Flory exponent v and the prefactor are related to the 
end-to-end distance of the polymer, i?oe = A^Nq. In the 
SS approximation, Eq. ([l} is solved simultaneously with 
Eq. ([3]), using T from Eq. ([3]) as a input in Eq. ([ij, and 
vice versa for s. 

For slightly smaller driving forces, 1 ^ / ^ Nq, the 
force is not sufficient to completely straighten the chain. 
Due to thermal fluctuations, a flower-shaped tail devel- 
ops (see Fig. [T]). In this stem-flower (SF) regime, the 
line density and and velocity of the monomers arc not 
constant in space. Therefore, one also has to solve the 
density an and the velocity vr near the tension front. As 
a result, one gets a system of equations. 



dR 
di 



= VR 



VqR , 

^ lyb ta.nh{b) 



cosh 6- 



R 



tanh (foa^/''-"V^ 



VR = Vq- 



tanh(6) 



VqR 



ln[cosli(fe)] 
6tanh(6) 



/tot - fjpVo 



(4) 
(5) 

(6) 
(7) 



that can be solved numerically for vq. Here, 6 is a (fixed) 
dimensionless parameter related to the spatial depen- 
dence of the velocity, and ensures global conservation of 
mass (see Appendix A). In the SF regime, ctq = 1, since 
the stem close to the pore is in a single-file configuration. 
The effective friction is given by Eq. ([2]). 

Finally, in the regime where the force insufhcient to 
straighten even a small part of the chain, / < , the 
chain adopts a trumpet-like shape. In this trumpet (TR) 
regime, the dynamics can be described by Eqs. dH)-®, 
with the velocity wo and density (Tq given by 



° &tanh(&) 
CTo = /tot - flpi'o 



/tot 

l-l/l/ 



l/u 



(8) 
(9) 



C. Pore friction 

The time evolution of the tension front (i? or A'') gives 
the contribution of the cis side subchain to the friction 




FIG. 2: (Color online) Comparison of waiting times w in both 
2D and 3D for MD and the BDTP model. The agreement of 
the BDTP model with MD simulations is excellent, and re- 
veals the two stages of translocation: the tension propagation 
stage of increasing w{S) and the tail retraction stage charac- 
terized by decreasing w{s). The parameters used were the 
same for both MD and BDTP (iVo = 128^^/ = 5, ksT = 1.2, 
= 0.7). The 3D MD resuhs are from [3g|. 



F. To complete the BDTP model, we still need to deter- 
mine the pore friction rjp. In general, rjp is a complicated 
function of the pore geometry, but here we restrict our 
study to the geometries used in our benchmark MD sim- 
ulations. In order to fix r/p, we examine the waiting time 
per monomer w(s), defined as the time that the individ- 
ual monomer spends inside the pore. With / sufficiently 
large, w oc F//. For small s, the friction F is mostly de- 
termined by f}p. Therefore, by comparing the w(s) of the 
BDTP model with MD simulations of Refs. [11, [2^, IH 
for the first few monomers, we have measured rj^ « 5 and 
7y2D 4 Jqj. ^Yic respective pore geometries. It should be 
noted that rjp is fitted only once, as opposed to being 
done separately for each combination of /, 77, etc. Thus, 
rjp is not a freely adjustable parameter. 



III. RESULTS AND DISCUSSION 

First, we validate the BDTP model through quanti- 
tative comparisons with MD simulations. In Fig. [5J we 
compare the waiting time if(s), which is the most im- 
portant and sensitive measure of the translocation dy- 
namics. As is shown, the match between BDTP and MD 
is almost exact. We stress that this agreement tells that 
the translocation dynamics is reproduced correctly at the 
most fundamental level and that such an agreement is a 
vital requirement for any correct theoretical model. The 
comparison also reveals an extremely lucid picture of the 
translocation process: first, as tension propagates along 
the chain, the effective friction increases and transloca- 
tion slows down. In the second stage, the number of 
dragged monomers is reduced as the tail retracts, and 
translocation speeds up. 

Next, we compare the exponents a obtained from the 
BDTP model with the corresponding numerical values 
from MD simulations. The parameter range has been 
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FIG. 3: (Color online) The effective exponent a(iVo) = 
in 2D (circles) and 3D (squares) as a function of A'o from 
the BDTP model, showing the extremely slow approach to 
the asymptotic limit a = \ ^r v . Most of experimental and 
simulation studies in the literature involve chain lengths of 
Afo < 10^ (shaded region), being clearly in the finite chain 
length regime. The inset shows the raw data t{No). Model 
parameters are the same as in Fig. (2] 

TABLE I: Values of a (r ~ A^o") from the BDTP model as 
compared to the corresponding values from MD simulations. 
a (BDTP) a (MD) Dimension and parameter values 

2D, T = 1.2, Ref. [18] 

1.51 ± 0.02 1.50 ± 0.01 / = 5.0, 7 = 0.7, 20 < % < 70 
1.71 ± 0.02 1.69 ± 0.04 / = 5.0, 7 = 0.7, 500 < iVo < 800 

1.52 ± 0.02 1.50 ± 0.02 / = 2.4, 7 = 0.7, 20 < No < 70 
1.71 ± 0.02 1.65 ± 0.04 / = 2.4, 7 = 0.7, 500 < A^o < 800 
1.66 ± 0.02 1.64 ± 0.01 / = 5.0, 7 = 3.0, 20 < A^o < 70 
1.71 ± 0.02 1.67 ± 0.03 / = 5.0, 7 = 3.0, 500 < A^p < 800 

3D, T = 1.2, Ref. [26] 
1.59 ± 0.02 1.58 ± 0.03 / = 0.5, 7 = 0.7, 16 < A^o < 128 
1.35 ± 0.02 1.37 ± 0.05 / = 5.0, 7 = 0.7, 16 < A^o < 256 
1.34 ± 0.02 1.37 ± 0.02 / = 10.0, 7 = 0.7, 16 < A^p < 256 

3D, T = 1.2, Ref. [25] 
1.41 ± 0.01 1.42 ± 0.01 / = 5.0, 7 = 0.7, 40 < A^p < 800 
1.39 ± 0.02 1.41 ± 0.01 / = 5.0, 7 = 0.7, 64 < A^p < 256 

3D, r= 1.0, Ref. [21J 
1.46 ± 0.02 1.47 ± 0.05 / = 3.0, 7 = 11.7, 70 < A^p < 200 
1.49 ± 0.02 1.50 ± 0.01 / = 30.0, 7 = 11.7, 200 < A^o < 800 



chosen to span the TR, SF and SS regimes, and to cover 
both short and long chain regimes in both 2D and 3D. 
The results are shown in Table I. The diversity of a in 
these regimes is evident, yet in all of them, the BDTP 
model is accurate to three significant numbers within the 
margin of error. This clearly shows, that while the values 
of a depend on several parameters, they all share a com- 
mon physical basis: non-equilibrium tension propagation 
on the cis side subchain. Then why is the exponent a not 
universal? The answer lies in the chain length regimes 
studied both in experiments and simulations. Typically, 
-^0 ^ 10'^ • However, in this regime a is not indepen- 
dent of the chain length! As shown in Fig. [3l a retains 
a fairly strong dependence on A^o up to A'o ~ 10^ in 2D 
and A^o ~ 10^ in 3D. Therefore, the observed scatter in a 
is a finite chain length effect, a fact that has been mostly 



ignored in the literature. 

Two additional remarks about the results of Table I are 
in order. First, the effect of the entropic term in Eq. ([1]) 
on a is extremely small. To show this, we solved Eq. ([T]) 
also without the entropic term. The results match ex- 
actly with those given in Table I, except for the low force 
case / = 0.5, r = 1.2 in 3D, where, without the entropic 
term, a = 1.56 instead of 1.59. Second, regarding / and 
Ao, the BDTP model gives two general trends for a: i) 
for a fixed /, a increases with A^o (as shown in Fig. [3]) 
and ii) for a fixed Nq, a decreases with /, as shown in 
Table I for 2D and 3D (Refs. [H,!!^). For f/ksT < 1, 
this trend is consistent with the value of a in absence of 
/ 0. For f/ksT » 1, a is almost independent of /. 
Therefore, the increase of a in the last two lines of Table 
I (Ref. [2l|) is in fact due to increase in Nq, not in /. 

Finally, we have estimated the asymptotic value of 
a in the SF regime by solving the BDTP model up to 
Ao = 10^° (Fig. E]). In 2D, the numerical estimate 
is « 1.750 ± 0.001, for A^o ^ 10^ and, in 3D, 

aj^ « 1.590 ± 0.002, for A^o > 10^. In both cases, we 
recover the value 1 -I- ^ as predicted in Ref. [3] with 
the constant-flux TP theory, and recently also using a 
different approach [l^. However, the value is differ- 
ent from Sakaue's orig inal prediction of ^+[+^'^' ^M^- 
The reason for the different predictions is that in both 
Refs. [illii, and in our model, the number of monomers 
is globally conserved, whereas in Refs. [TT| - [l3j the con- 
servation is guaranteed only locally in the neighborhood 
of a; = —R. Therefore, asymptotically, a^o ~ 1 + in 
both 2D and 3D, also in agreement with the prediction 
of Ref. [8i]. 



IV. CONCLUSIONS 

To summarize, we have introduced a new theoretical 
model of driven polymer translocation that has only two 
degrees of freedom and no free parameters. The model 
gives near-exact agreement with high-accuracy molecular 
dynamics simulations in a wide range of parameters. Our 
study shows that the dynamics of driven translocation is 
dictated by non-equilibrium tension propagation on the 
cis side subchain. The model also reveals that the major- 
ity of experiments and simulations in the literature are 
performed in the regime, where finite chain length effects 
have significant impact on the translocation dynamics. 
Although mostly overlooked in the literature, this is an 
important observation, since the finite chain length ef- 
fects persists for chain lengths of at least several tens of 
thousands of monomers. Therefore, in most studies of 
polymer translocation, finite chain length effects cannot 
be neglected. This fact is also vital for the theoretical 
study of driven polymer translocation. 



5 



Acknowledgments 

This work has been supported in part by the Academy 
of Finland through its COMP Center of Excellence and 
Transpoly Consortium grant. TI acknowledges the finan- 
cial support of FICS and TES. The authors also wish to 
thank CSC, the Finnish IT center for science, for alloca- 
tion of computer resources. 



Appendix A: Derivation of the tension propagation 
equations 

For convenience, we use dimensionless units denoted by 
the tilde symbol as X = X/Xu, with the unit of length 
flu = a, force /„ = kgT/a, time t„ = rjo? /ksT, veloc- 
ity Vu = a/tu and friction rju = t], where 77 is the solvent 
friction per monomer. The tension propagation (TP) for- 
malism is derived for A^o 00 in Refs. [Il| - [l4| . Here, we 
derive the formalism for finite chain lengths by including 
the explicit pore-polymer interactions through the pore 
friction fjp, and the spatially dependent velocity profile 
for the polymer chain. These will be discussed below. 

To begin, we note that as the driving force at the pore 
is applied, the chain starts to move in stages, as tension 
propagates along the backbone. At time t, when s(t) 
monomers have translocated, N{t) — s{t) monomers on 
the cis side are under tension and moving towards the 
pore at velocity v(x,t) and, the remaining Nq — N{t) 
monomers are at rest. The moving and unmoving do- 
mains are separated by a tension front at i = —R. 
The subchain between — i? < i < is deformed un- 
der the tension and can be regarded as a self-avoiding 

The 



walk following the Pincus blob description [31 



subchain adopts a configuration of increasing blob radii 
^(i), with the blob closest to the front having radius 
= ^{—R + £,r)- At length scales shorter than the blob 
size £^{x) = the chain behaves as if undisturbed 

by the driving force, scaling as ^ = 9'^ 1 where g is the 
number of monomers inside the blob and v is the Flory 
exponent. This gives the relation a = g/£^ = for 
the monomer line density a. By definition, tr satisfies 
I-R{i} ^(^'^ i)dx' = N{i) — s{i). To solve i{x), we require 
the local force balance between the driving force and the 
drag force at x: 



v{x\ imx', t)Yi''-Ui' = f{i, t) = t)-\ 

R{i) 

(Al) 

In addition, there is a balance between the driving force 
at the pore entrance, /o = /(O) and the total drag force 
of the cis side subchain. This global force balance can be 
enforced by substituting i = to Eq. (jAip : 



The time evolution of the tension front R obeys the 
equation of conservation of monomers. 



dR{i) 
di 



Mi) 



dN{i) 
di 



(A3) 



Here we employ the short-hand notation aR{t) = a{—R+ 
^R,i) and VR{i) = v{—R + £,R,i). Furthermore, since the 
monomers outside the tension front arc on the average 
immobile, the location of the tension front R is given 
by the equilibrium end-to-end distance of the subchain 
consisting of the first N monomers [Tl| : 



(A4) 



Here, A^, is a model-dependent prefactor for a chain 
with one end tethered to a wall. For the ideal chain, 
Al, = 1 and, for the self-avoiding chain we have mea- 
sured « 1.16 ±0.05 in 2D and 1.15 ±0.03 in 3D from 
MD simulations using the Kremer-Grest model [3I] with 
typical parameter values for the chain. Details of the 
simulation are explained in Appendix B. 

To couple the TP model with the effective friction of 
the Brownian dynamics equation of s, we write down 
the conservation of monomers at the pore entrance, 

^ = (To(i)wo(i), which also defines f(i) = aoi)t,{i) ' 
where <7o{i) = (7{0,i) and vo{i) = v{0,i). Finally, we 
define the explicit relationship between the total driving 
force /tot and the force at the pore entrance, fo- For fi- 
nite A^Oi we need to take into account the pore-polymer 
interactions by introducing the pore friction coefficient 
Vp = Vp/V- This defines the relationship between /tot 
and fo as /tot - VpVQii) = fo. Obviously, for an ideally 
frictionless pore, rip — and /tot = fo- For finite A^o and 
77p, the overall effect of fjp is to make the translocation 
time T less sensitive to No, i.e., decrease a. 

Strong stretching (SS) regime. Let us first consider 
the case of a strong driving force so that the moving part 
of the chain is almost completely straight. This strong 
stretching regime is realized when fo > [U . In the 
SS regime, the line density of monomers in the moving 
domain is constant d'{x) ~ a* ~ I and correspondingly 
S,{x) = ^* « 1. As a first approximation, one may assume 
that the segment length a remains unchanged. Although 
a more detailed treatment is possible [Tsj , this is a reason- 
able approximation in the relevant range of forces. From 
these assumptions, one immediately obtains a constant 
velocity profile for the monomers in the moving domain: 
v{x,i) = vo{i)Q{x + R), where 6(i) is the Heaviside 
step function. In particular, vo{t) = VR{t) Vi. In the SS 
regime, one obtains the TP law 



dN 
di 



fu 



r(<)(i-z/A,iv--i)' 



(A5) 



Rii) 



v{x ,t)a{x ,i)dx = fo 



(A2) 



with r{t) = N{t) - s{t) + f)p. Note that Eq^^ (|A5|) is a 
generalization of the TP law of Refs. [ij, [T^ to pore- 
driven translocation with an explicit pore friction rjp. 
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Trumpet (TR) regime. In the opposite limit of ex- 
tremely weak driving force (/o i$ ^o'^)i the whole chain 
can be considered to be in equilibrium, and one recov- 
ers the Rouse-type friction of the subchain. However, 
such small forces arc rarely realized in simulations or ex- 
periments. Indeed, if the force is only slightly larger, 
■^o''^ ^ /o ^ 1 ■ the chain adopts a configuration resem- 
bling a trumpet, where the blob radius f_increases as one 
moves further away from the pore [ll|, . This leads to 
(To(?) < S-R{i) Vf. In Refs. j7l| - [l3j . it is assumed that the 
velocity profile is constant in the TR regime, similarly 
to the SS regime. This assumption, however, leads to a 
contradiction. By integrating Eq. (|A3|) over the whole 
translocation process from i = to i = f, it follows that 
s{f) < Nq. This is an obvious contradiction, a point that 
was also raised in Ref . • The only way to remove the 
contradiction is to relax the constant-velocity assump- 
tion so that vo{i) > vji{i). It should be noted that the 
constant-flux approximation of Ref. jisl is also not valid 
in the short chain regime, where dR/dt is of the order of 
ds/dt. To solve this problem, we have studied the veloc- 
ity profile using MD simulations, and find that at least 
for iVo < 10'^, the velocity profile is to good approxima- 

~/~ ~\ ~ tanh\b(x/ R+l]\ [7—1 

tion given by v[x,t) = voit) jj^^b) — - 135|J. Here b is 

a dimensionless parameter that controls the sharpness of 
the profile and is fixed by enforcing global conservation of 
monomers, i.e., requiring that s(f) = Nq and R{f) — 0. 
Because the approximate profile for v{x, t) is not exact, b 
has a weak dependence on chain length iVo and is solved 
for each iVo by numerical iteration until s(f) = iVo and 
R{t) = is satisfied. 

In the TR regime, the line density a of the monomers 
is not fixed, but is determined by Eq. (|Aip . The line 
density near the pore is determined by the blob radius 
^(0) = fo^- To calculate the line density at the bound- 



ary, one has to solve Eq. (|Aip for x G 



-R,-R + ^f 



The resulting implicit equation for is solved numeri- 
cally: 



i^6tanh(6) 



In 



cosh 

[ Rit) , 



(A6) 



The velocity v^it) can be solved from Eq. (|Aip with x = 
0, which reduces into another implicit equation. 



~ /rNp/rNln[cosh(6)] 
VQ{t)R{t)- 



/tot - j)pUo(i) 



6tanh(6) 

The velocity near the boundary is given by 

tanh (biR/R^ 



ijR{t)=i{-R + ^n,t) = io{t)- 



tanh(6) 



(A7) 



(A8) 



The equation of motion for the tension front can be solved 
from Eqs. (|A3l) and (jA4|) that give 



dR{i) 
di 



iRii) 



-AZ^I''dR{t)-^R{tf'''-^ 



(A9) 



Finally, the time evolution of the tension front is given 
by numerically solving Eqs. (|A6p - (|A9p with the initial 
condition R{tq) « f^^ 



Stem-flower (SF) regime. In the intermediate regime, 
1 ^ /o ^ ^0 J the chain assumes a shape consisting of a 
fully elongated stem between the pore and x = —r, fol- 
lowed by a trumpet-shaped flower for —R < x < —f. 
The analysis of the stem is similar to the SS regime, 
and the flower in turn follows the TR regime calculation. 
The two parts are connected via the boundary condition 
/(— f) = 1, which, together with Eq. ()Aip evaluated at 
— f, can be used to eliminate f. giving the equation for 
Mi)- 



,~, ~ ,~, ln[cosh(fe)l 



ftot-VpVa{i) +1^-1. (AlO) 



The blob radius and the velocity at the boundary are 
given by Eqs. (|A6|) and (|A8|) . respectively, and the time 
evolution of the front by Eq. l(M\i . Note that Eqs. ([Tf)) 
and (|A10[) ensure a smooth cross-over between the TR 
and SF regimes at /o = 1 and that the SF regime equa- 
tions approach the SS regime Eq. (jA5|) when / ^ 1 
(as f — >■ i?, — > 1 and b — > oo). In practice, we 
solve Eqs. ((X3|) and (|M1) - (|A10|) . choosing Eq. (|T7| over 
Eq. (|A10p if /o < 1, and vice versa. 



Appendix B: Molecular dynamics simulations 

The details of the molecular dynamics simulations that 
wc have used for benchmarking the BDTP model arc ex- 
plained in this Section. In the MD simulations, the poly- 
mer chain is modeled as Lennard- Jones particles inter- 
connected by nonlinear FENE springs. Excluded volume 
interaction between monomers is given by the short-range 
repulsive Lennard- Jones potential: 



f/Lj(r 




e for r < 2^/^ 
for r > 2i/V 



(Bl) 



Here, r is the distance between monomers, a is the diam- 
eter of the monomer and e is the depth of the potential 
well. Consecutive monomers are also connected by FENE 
springs with 

;7FENE(r) = ^hRl ln(l - r^R^, (B2) 

where k is the FENE spring constant and Ro is the maxi- 
mum allowed separation between consecutive monomers. 
For the chain, we use the parameters e = 1, fc = 15 and 
Rq ~ 2. The main part of the wall is constructed us- 
ing a repulsive external potential of the Lennard- Jones 



4e 



(f)"-(fr 



e for x < 2^/^(7 and 



form C/cxt 

otherwise. Here x is the coordinate in the direction 
perpendicular to the wall, with x < signifying the cis 
side and a; > the trans side. The neighborhood of the 



7 



pore is constructed of immobile Lennard- Jones beads of 
size a. All monomer-pore particle pairs have the same 
short-range repulsive LJ-interaction as described above. 
We have verified that using the simple external potential 
C/oxt for the wall gives the same results (within statistical 
error) as using a wall made of monomers in fixed lattice 
sites, at least as long as the interaction between the wall 
and the polymer is purely repulsive. 

Similarly to most of the molecular dynamics simula- 
tions in the literature, we take the surrounding solvent 
into account through frictional and random forces. Thus, 
each monomer is described by the Langevin equation of 
motion 

mi=j = -V([/l.i + t^FENE + Uc^t) - ?7Vj + Q-, (B3) 
where m is the monomer mass, 77 is the friction coef- 



ficient, Vi is the monomer velocity, — Vf/ext = / is the 
external force in the pore and C,i is the random force with 
the correlations {Ci{t)Cj{t')) = 2r]kBT6ijS{t — t'), where 
ksT is the thermal energy. Typically, we have used the 
parameter values m = 1, 77 = 0.7, ksT = 1.2. The equa- 
tions of motion are solved with the BBK algorithm [3^1 
with time step St = 0.005. Initially, the polymer chain is 
placed on the cis side with the first monomer located at 
the pore entrance. Eq. ()B3|) is then solved numerically 
while keeping the first monomer fixed until an uncorre- 
latcd initial configuration is generated. After that, the 
whole chain is allowed to evolve according to Eq. (|B3[) un- 
til the chain escapes either to the cis or trans side. The 
latter is recorded as a successful translocation event. We 
average our MD results over at least 2000 such events. 
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